> #Riadioaktyviosios medziagos skilimo modeliavimas
> # Kintamuju reiksmes
>  #alpha - skilimo greiti aprasanti konstanta
>  #P - procentai, radioaktyviosios anglies likutis iskasenoje
>  #t - laikas(metais)
>  #m - mase (gramais)
>  #m0 - pradine mase (gramais)
> restart;
> DL := diff(m(t),t) = -alpha*m(t);
                         d                    
                  DL := --- m(t) = -alpha m(t)
                         dt                   
> Bendrasis_sprendinys := dsolve(DL,m(t));
        Bendrasis_sprendinys := m(t) = _C1 exp(-alpha t)
> #Radioaktyviosios anglies kiekis medziagoje  per 5730 metu sumazeja per pus, taigi galime apskaiciuoti skilimo greitiaprasancia konstanta (alpha)
> t0 := 5730:
> F:=(ln(1/y)/t):
> Apskaiciuojam_alpha := evalf(subs(y=2,t=t0,F),5);
               Apskaiciuojam_alpha := -0.00012097
> Ieskoma_lygtis:=subs(alpha = Apskaiciuojam_alpha,c=c0,Bendrasis_sprendinys);
                                   10                        
   Ieskoma_lygtis := 4.600000000 10   = _C1 exp(0.00012097 t)
> Laikas:=(ln(P)/Apskaiciuojam_alpha);
                  Laikas := -8266.512358 ln(P)
> m0 := 5.0*10^10:
> m := 4.6*10^10:
> Kiek_procentu_pasilike_medziagoje := simplify(m/m0);
       Kiek_procentu_pasilike_medziagoje := 0.9200000000
> Lygtis := (ln(Kiek_procentu_pasilike_medziagoje)/(Apskaiciuojam_alpha)):
> Iskasenos_amzius_(metais) := evalf(Lygtis,3);
               Iskasenos_amzius_(metais) := 689.
> p:=(exp(Metai*Apskaiciuojam_alpha)*100):
> plot(p, Metai=0..30000, thickness = 1, color = blue);

> 
